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ABSTRACT 

Context. The role of magnetic fields in gamma- ray burst (GRB) flows remains debated. If strong enough, they can leave their signature on 
the initial phases of the afterglow by substantially changing the backreaction of the flow as a consequence of its interaction with the external 
medium. 

Aims. We aim to quantitatively understand the dynamical effect and observational signatures of magnetization of the GRB ejecta on the onset 
of the afterglow. 

Methods. We perform ultrahigh-resolution one-dimensional relativistic MHD simulations of the interaction of a radially expanding, magne- 
tized ejecta with the interstellar medium. The need of ultrahigh numerical resolution derives from the extreme jump conditions in the region 
of interaction between the ejecta and the circumburst medium. We study the complete evolution of an ultrarelativistic shell all the way to a the 
self-similar asymptotic phase. 

Results. Our simulations show that the complete evolution can be characterized in terms of two parameters, namely, the £ parameter introduced 
by Sari & Piran (1995) and the magnetization cr . We exploit this property by producing numerical models where the shell Lorentz factor is 
70 ~ tens and rescaling the results to arbitrarily large yo. We find that the reverse shock is typically very weak or absent for ejecta characterized 
by cro ~ 1- The onset of the forward shock emission is strongly affected by the magnetization. On the other hand, the magnetic energy of the 
shell is transferred to the external medium on a short timescale (~several times the duration of the burst). The later forward shock emission 
does not contain information for the initial magnetization of the flow. The asymptotic evolution of strongly magnetized shells, after they have 
suffered a substantial deceleration, resembles that of hydrodynamic shells, i.e., they fully enter in the Blandford-McKee self-similar regime. 
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1. Introduction 

It is generally believed that gamma ray bursts (GRBs) are the 
result of the energy release in an ultrarelativistic outflow. The 
mechanisms responsible for launching, accelerating the flow 
and powering the GRB emission are still not well understood. 
Two alternatives for the energy content of the flow have been 
explored. The GRB flow may be in itially dominated by thermal 
energ y density creating a fireball ( iGoodman Paczynski 
19861) or by the energy stored in magn etic fields giving rise to 



a Poynting-flux dominated flo w (PDF ; lUsovll 1 992t iThomp son 
19941: iMeszaros & Rees|[l997l) . 

In fireball models, magnetic fields are not dynamically im- 
portant at any stage of the flow evolution. On the other hand, 
models of MHD jet acceleration indicate that the conversion 
of Poynting flux to kinetic energ y is partial dMichell [l969t 



Li et all 1 19921: iBeskin et all 1 19981: Iviahakis & Konigll 12003 



Giannios & SpruiJ 20061) . As a result, the flow remains rather 
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strongly magnetized at large distance where it interacts with 
the circumburst medium. 

The interaction of the fast flow with the external medium 
likely powers the afterglow emission. The initial phases of this 
interaction depend, to a large extent, on the magnetization of 
the flow. Strong magnetic fields affe ct the shock condition s and 
the internal dynam i cs of the ejecta. iRees & Gunnl ( 119741) and 
Kennel & Coronitil ( 119841) have paved the way to understand- 
ing the effect of the magnetization on the shock conditions. 
They showed that, in the ideal MHD limit, shocks in mag- 
netically dominated flows cannot efficiently dissipate energy. 
This particular effect has been studied recently by examining 
the shock conditions of the (reverse) shock that propagates 
into the GRB flow and the resulting emission dFan et alj|2004c 
Zhang & Kobavashil2005 ). Ta king the internal evolut ion of the 
magnetized shell into account. iGiannios et alj (2008) have ar- 
gued that even moderately strong magnetic fields can suppress 
the reverse shock altogether possibly explaining the observed 
paucity of reverse shock signature in the early afterglow. 
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Here, we complete previous (semi-)analytical studies of the 
afterglow phase of magnetized ejecta with relativistic MHD 
simulations. We follow the deceleration of magnetized ejecta 
from the initial phases of the interaction up to the self sim- 
ilar regime where all the energy has been transferred to the 
shocked external medium. These simulations clarify what are 
the dynamical effects of magnetization of the GRB ejecta and 
their observational implications connected to the forward and 
reverse shock emission. Along the way, we will show a new 
set of scaling laws jj4.4l that enable us to extrapolate the results 
of numerical models with moderate values of the initial bulk 
Lorentz factor (~ 15) of the ejecta to equivalent models with 
much larger Lorentz factors >100. 



2. Ejecta-medium interaction 

At large distances from the central engine, there is a substantial 
interaction of the relativistic ejecta with the external medium. 
This interaction is believed to result in the afterglow emission. 
An important difference between fireballs and PDFs lies in the 
magnetization of the ejecta at the onset of the afterglow phase. 
In fireball models, the energy of the flow is dominated by the 
kinetic energy of baryons at large distance from the central en- 
gine. If the flow is launched Poynting-flux dominated, it is ex- 
pected to maintain a large fraction of its energy in the form of 



magnetic energy, the rest being in kinetic form (IMiche 



Li et al.lll992t Beskin et al.l ll998l: bre nkhahn & Sprui 



it 



1969; 



2002; 



Lyutikov & Blandfordll2003l) . Since (at least the initial phases) 
of the ejecta-external me dium interaction depend o n the mag- 
Fan et al 



n depend o n tf 

netization of the ejecta dKennel & Coronitil Il984t |Fi 

et al.ll2007T> . it 



20041 IZhang & Kobayashill2005l iGene" 



is pos- 
sible to discriminate among fireballs and PDF models from af- 
terglow observations. As we show in this work, early afterglow 
observations are particularly promising in this respect. 

2.1. Previous Studies 

The deceleration of non- magnetized ejecta has been well stud- 



ied w ith both analytical dRees & Meszarosll 1992c I Sari & Piran 



995l) and with numer ical approaches in one di mension (ID) 
(Kobayashi et al .1 19991) and two dimensions (2D) ( IGranot et al 



2001tlMeliani et alj2007l) . The 2D studies are importa nt to fol 



low th e late-time lateral spreading of collimated ejecta (Rhoads 
1999b . On the other hand, the initial phases of the deceleration 
of the ejecta, in which we are interested here, are not affected 
by 2D effects and can be studied assuming spherically symmet- 
ric flow. 

Sari & Piranl ( 1995 ) considered the case of non-magnetized 



ejecta assuming a cold shell with (isotropic equivalent) kinetic 
energy E, Lorentz factor yo and width Ao, which moves against 
external medium with density p e . This interaction leads to a 
pair of shocks: one that propagates in the external medium 
(forward shock) and one that slows down the ejecta (reverse 
shock). The strength of the reverse shock depends on the ratio 
of the densities of the shell and the external medium and on 
the bulk Lorentz factor of the flow. It can be shown that the 



strength of the reverse shock can be conveniently parametrized 
by 



I 1 



(i) 



where / = (3£ , /4mi e »JpC 2 ) 1/3 is the Sedov length, n e is the num- 
ber density of electrons in the external medium and m p is the 
proton mass. E is the total energy (kinetic in this case) of the 
ejecta. In the limit where £ » 1 the reverse shock is Newtonian 
and the shell is said to be "thin". The ejecta do not decelerate 
much by the time the reverse shock crosses them. If £ <K 1, we 
find ourselves in the "thick shell" case, and the reverse shock is 



relativistic and slows down the ejecta appreciably (Sari & Piran 
1995J). . For typical parameters of GRB flows £ is of order of 
unity with 0.1 < f< several. 

The dynamics of the deceleration of strongly magnetized 
ejecta have not been studied in the same detail. In addition to 
the £ parameter, the ejecta are characterized by the magneti- 
zation er g defined as the ratio of m agnetic-to-kinetic energy in 
the flow. iKennel & Coronitil d 19841) solved for the ideal MHD 
shock conditions for arbitrarily magnetized ejecta with a domi- 
nant toroidal field, and showed t hat the dissipation by the shoc k 
gets weaker as <tq increases (e.g lLvutikov & Blandfordll2003l) . 
They applied their a nalysis to the standing shoc k of pulsar 
winds. More recently IZhang & Kobayashil (12005b focused on 
the effect of magnetization in the context of GRB afterglows. 
They ignored the internal evolution of the shell prior to the 
interaction with the external medium and studied the reverse 
shock crossing phase (provided that there is a reverse shock 
forming). They found distinct features in the early time light 
curves because of the magnetization. The shocks from the in- 
teraction of the GRB ejecta with the external medium propa- 
gate forwards and the shock conditions depend on the di stance 
from the central engine. The Zhan g & Ko bayashil d2005l) anal- 
ysis has been criticized by iLyutikovl (120051) f or the assump 



tion on the distance where the ejecta decelerate. iGiannios et al 



(120081) took into account the internal evolution of the ejecta 
and derived the an analytic condition for existence of a reverse 
shock depending on £ and cr in a parameter space relevant for 
GRB flowsQ 

After the reverse shock (if there is one) reaches the back 
part of the ejecta, there is a transient phase of interaction where 
rarefaction waves cross the shocked ejecta and shocked exter- 
nal medium. Gradually most of the energy is passed in the 
shocked external medium and the whole structure relaxes to 



the self-similar blast wave described in Blandford & McKee 



(119761) . From this point on, the evolution of the blast wave de- 
pends only on the total energy E and the density of the external 



1 At the time of the referring of this work, iMizuno et alj (2008) 
have also published a work addressing the problem of the deceleration 
of arbitrary magnetized ejecta into an unmagnetized medium, and dis- 
cuss its implications for GRBs and active galactic nuclei. However, in 
spite of the undoubtable academic interest of their studies, the con- 
ditions set by these authors (particularly, the use of planar symmetry, 
and the very small density contrast between the magnetized shell and 
the external medium ~ 100) are far from those met in GRB afterglows, 
specially during the early afterglow propagation. 



Mimica, Giannios & Aloy: Deceleration of arbitrarily magnetized GRB ejecta 



3 



medium « e and not on ctq. After the self-similar evolution has 
been reached, nothing can be inferred about the initial magne- 
tization of the flow. 

However, none of these studies have addressed two impor- 
tant questions. First, there is a question in which stage of the 
interaction a reverse shock forms (if it forms at all). At short 
distance from the central engine the magnetic pressure of the 
shell is high enough that the shell rarefies upon interacting with 
the external medium. This rarefaction may turn into a (reverse) 
shock at larger distance where the magnetic pressure in the 
shell drops. Second, although it is clear that (ignoring radia- 
tive losses) the total energy initially in the shell is passed onto 
the external medium at a distance of the order of the Sedov 
length which is independent of the magnetization of the flow 
dLvutikovll2005l) . the details of how exactly this happens de- 
pend on the magnetization. These two aspects are closely con- 
nected to the energy that is dissipated in the (forward and re- 
verse) shocks as function of distance and, consequently, to the 
afterglow emission from particles accelerated in these shocks. 
We address this issue here by studying the full dynamical in- 
teraction from the initial stages all the way to the establish- 
ment of the self-similar evolution. To this end we perform ultra- 
high resolution one-dimensional relativistic MHD simulations 
of shell-medium interaction. 



3. The model for the ejecta 

We focus on the GRB flow at a distance where there is sub- 
stantial interaction with the external medium. This interaction 
likely takes place well after the acceleration, collimation and 
prompt emission phases are over. After the internal dissipation 
phase (believed to power the y-ray emission) finishes, the flow 
expands radially and cools down. The expansion also leads to 
a dominant toroidal component for the magnetic field. At the, 
so-called, Alfven point, the poloidal B p and toroidal field 
components are expected to be of similar magnitude. Further 
out, the flux freezing condition results in B p oc 1 /r 2 while the 
induction equation predicts slower decline for oc 1/r. The 
same scalings hold if the initial B-field were random resulting 
in Z?^ » Z? p at large distance from the central engine. 

We consider radially moving, cold shell of ejecta of width 
Ao, total (kinetic and magnetic) energy E that coasts with 
a bulk Lorentz factor yo. The magnetic content of the flow 
is parametrized with the magnetization parameter <xo which 
stands for the magnetic-to-kinetic energy ratio in the shell. The 
flow is assumed to move with super-fast magnetoso nic speeds 



(i.e., yi > 1 + cp; for stu dies of the opposite limit see lLyutikov 



200dlGenetetall2007h . For the simulations presented bellow, 
the shell is located at an initial distance ro from the central en- 
gine. The choice of ro is important since ro should be small 
enough not to affect the subsequent interaction of the ejecta 
with the external medium, ro must be set smaller than any of 
the characteristic radii that appear when considering the de- 
celeration of magnetized ejecta. These radii are the 'contact' 
radius and the 'reverse shock crossing' radius to be defined in 
the next section. 



3.1. Characteristic distances 

In the super-fast magnetosonic flow under consideration, the 
various parts along the radial direction have dropped out of 
MHD contact during the acceleration phase. It can be shown 
that for a cold flow that coasts with constant speed with dom- 
inant toroidal field the magnetization remains constant. The 
time it takes for a fast MHD wave to cross the width of the 
shell is therefore fixed. The expansion timescale f exp = r/yoc 
is initially much shorter than that of MHD waves but in- 
creases linearly with distance from the central engine. At the 
so-called 'contact' radius r c MHD waves cross the width of 
the shell on a timesca le comparable to the expansion timescale 
dGiannios et al.ll2008l) 



AoTo 



(2) 



After contact is established, the shell is no longer in pure ballis- 
tic motion and internal evolution because of MHD forces can 
no longer be ignored. On the other hand in non-magnetized 
ejecta the sound speed drops fast with distance because of adi- 
abatic expansion, and the motion is is not affected by the pres- 
sure of the shell. 

A second important radius is the radius where the reverse 
sh ock reaches the rear part of the ejecta. This radius is derived 
by lZhang & Kobayashil (120051) by solving the ideal M HD shock 



condi tions for arbitrarily magnetized ejecta (see also lFan et al 



2004 for the case of mildly magnetized ejecta). Their analysis 
describes the reverse shock crossing phase provided that there 
is a reverse shock forming. The re verse shock crossing radius 
can be approximately expressed as (IGiannios et al 1 l2008l) 



r rs ^X /4 /VlT^o"- 



(3) 



The IZhang & Kobayashil analysis does not take into account 
the internal evolution of the shell. It is thus accurate when such 
evolution is not significant, i.e. when r rs < r c . 

The initial distance ro where the shell is set up must be 
ro <K min[r c ,r rs ] so that the simulation starts early enough to 
follow both any rarefaction waves within the shell, and shock 
waves result from the interaction with the external medium. 

3.2. Characteristic quantities 

In this paper we frequently use the following definition of the 
normalized time of observation for a model with parameters yo, 
Ao and ro: 

fobs := Aq 1 [ct - r] , (4) 

where f b s is the time of observation of a signal sent from radius 
r at time t in the GRB frame or laboratory frame, normalized 
to the light crossing time of the initial width of the shell Ao/c. 
As we will show in Sec. 14.4.31 this definition of f 0DS enables us 
to compare properties (in the observer frame) of shells with the 
same f independent of their initial Lorentz factor. 

We also often base our arguments on the relative Lorentz 
factor y re i between two parts of the fluid separated by the shock 
front. For ultrarelativistic flows we use 

1 (7a j_ Jb \ 

y re i:=- — + — , (5) 
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Table 1. Parameters of the numerical models 



O"0 





1 


3 


thin shell (£ = 1.1) 


V 


V 




thick shell - 0.5) 


V 


V 




continuous flow 




V 


V 



where y a and yj are the Lorentz factors of the fluid ahead and 
behind of the shock, respectively. We point out that y re i depends 
only on the ratio y a /y&. 

3.3. Numerical models 

Although the problem is characterized by several parameters E, 
n e , Aq, yo and o"o, it turns out that for the systematic study of the 
shell-medium interaction we need to focus on the combination 
of the first four ones parametrized by 



^(£/« e ) i/6 /Ai /2 r r 



(6) 



and ctq. We demonstrate and quantify this statement in the next 
Section. In order to simplify the analysis, we restrict ourselves 
to the case in which the external medium density is uniform 
(p e 3 x l(T 4 po(i + o"o); po being the initial shell density, and 
leave the study of stratified external media for a future work. 

As in the case of unmagnetized ejecta, we make use of 
the Sari-Piran parameter £ in order to partly characterize the 
strength of the reverse shock. Certainly, in the magnetized 
case, the shock strength is not uniquely set by Instead, 
an additional parameter, <x , needs to be introduced to fully 
describe the reverse shock strength of arbitrarily magnetized 
flows. Thus, one deals with a £ - <x plane in exploring dif- 
ferent cases for the initial phases of shell-external medium in- 
teraction. Here, we explore the £ ~ 1 regime that is relevant 
for typical GRB parameters. Numerical reasons limit us to the 
< o"o < 3 range for the magnetization parameter. 

Our model runs are summarized in Tab. (T| . The £ = 1.1 
runs (thin shells) are characterized by E = 3.33 x 10 53 erg, A = 
10 15 cm, ro = 15, « e = 10 cm -3 . The £ = 0.5 models (thick 
shells) have ten times larger total energy E and width of the 
shell Aq while yo and n e remain fixed. The "continuous flow" 
model (to be discussed in more detail in Sect. 4.3) describes 
a flow of constant total (kinetic+Poynting) luminosity of L = 
10 49 erg/sec that moves with yo = 15 and collides with external 
medium of number density n e = 10 crrT 3 . In all models ro = 
5 x 10 16 cm. 

One may notice that the model runs are characterized by 
unrealistically low Lorentz factor yo =15 and wide shells Ao ~ 
10 15 cm with respect to what is expected from a GRB flow (i.e. 
yo ~ 100, Ao < 3 x 10 12 cm). This choice of parameters is made 
so that the problem is reliably resolved with our RMHD code. 
While runs with yo > 100 in combination with extreme density 
and magnetic field jumps at the edge of the ejecta shell are not 
feasible at this stage, we propose a method to extrapolate the 
results of the yo = 15 simulation by appropriately rescaling 
of the initial conditions. Furthermore, we have run simulations 
with yo = 10,20 where we demonstrate the accuracy of the 
rescaling procedure (see Section 4.4). 




r [10 17 cm] 

Fig. 1. Snapshot of the thin magnetized shell evolution taken 
after the RS has formed and before it has crossed the shell. Full 
and dashed black lines show the logarithms of the rest-mass 
density (normalized to the initial shell density po) and of the 
pressure (normalized to poc 2 ). The red line shows the logarithm 
of the magnetization <x, while the blue line shows the fluid 
Lorentz factor y in the linear scale. All quantities are shown 
as a function of radius r. Positions of the forward shock (FS), 
contact discontinuity (CD), reverse shock (RS), and left- (Ri) 
and right-going (R2) rarefactions have been indicated. There 
exist another rarefaction (moving backwards in the external 
medium) and another contact discontinuity that form at the rear 
edge of the shell, which are not shown here. Both structures are 
located to the right of Ro- The unshocked external medium is 
located in front of the FS and its density p e w 6 x 10~ 4 po is 
smaller than the minimum density shown in the plot. The rar- 
efaction R2, the CD and the FS display a qualitatively similar 
profile in the non-magnetized case. The rarefaction Ri and the 
late steepening of the conditions at its tail, resulting into the 
formation of RS, are specific of magnetized ejecta. 

4. Results 

To derive the results presented in this Section, we solve the 
equations of RMHD in ID spherical geometry with mag- 
netic field perpendicular to the direction of propagation of the 
fluid, i.e., with a purely toroidal magnetic field. The system 
of RMHD equations, and the numerical tests we have made 
to choose the appropriate numerical resolution for our experi- 
ments (between 10 4 and 6 x 10 4 cells to resolve the initial radial 
width of the ejecta) are shown in the Appendix [A] Finally, all 
our models have been run until the bulk Lorentz factor behind 
the forward shock has dropped to y ~ 2 - 3. By that time the 
shell has suffered a substantial deceleration and fully entered in 
the Blandford-MacKee self-similar regime. 



4.1. Non-magnetized shells 

The 0" O = models sh ow the well known features expected 
from analytical results ( I Sari & Piranl 1 1995b and simulations 
( Kobayashi et al 1 fl999n . The thin shell (£ = 1.1) model is char- 



Mimica, Giannios & Aloy: Deceleration of arbitrarily magnetized GRB ejecta 



5 



3-1 



3 

M 
O 



2-1 



M 
O 




r [cm] 



Fig. 2. Energy in different components (normalized to the ini- 
tial total energy in the shell) as a function of the radius of the 
FS. Upper and lower panels show the thin = 1.1) and thick 
= 0.5) shell models, respectively. Black and blue lines cor- 
respond to un-magnetized (<x = 0) and magnetized (cr = 1) 
models, respectively. Thick full, dashed and dotted lines re- 
spectively show the kinetic, thermal and magnetic energy of the 
shell. The thin full line shows the total energy in the shocked 
external medium. Thick vertical lines denote the radius of the 
FS at the moment when the RS has crossed the shell. 



obs 



Fig. 3. Relative Lorentz factor at the FS (upper panel) and RS 
(lower panel) as a function of the normalized observer time f b s 
(EqH|). Full and dashed lines correspond to the thin and thick 
shell models, respectively. Black and blue colors denote non- 
magnetized and magnetized models, respectively. The vertical 
lines in the lower panels denote the time of observation when 
the RS crosses the shell. In the upper panel the vertical lines 
denote the time of observation when the rarefaction R2, that 
originates from the rear edge of the shell and interacts with the 
RS, catches up with the FS. 



acterized by a Newtonian-to-mildly-relativistic reverse shock. 
The reverse shock crosses the shell at a distance r ls = 3 x 10 17 
cm which agrees within ~ 10% with the analytical estimate 
from Eq. ©. At this distance, the relative Lorentz factor of the 
shocked ejecta with respect to the unshocked shell is y re \ 
1.18. The thick shell (£ = 0.5) model finds itself closer to 
the "relativistic reverse shock" regime with r rs = 9.5 x 10 17 
cm (within ~ 10% of the initial estimate) and y re i 1.60. In 
both runs, after the reverse shock crosses the shell, there is a 
rarefaction that starts at the rear part of the shell and prop- 
agates forwards (a similar rarefaction happens in the magne- 
tized case, which we label R2 in Fig.Q]i. The rarefaction crosses 
the contact discontinuity, generated in the leading radial edge 
of the shell (an equivalent contact discontinuity arises in the 
magnetized case, see "CD" in Fig. [TJ, and reaches the forward 
shock when the shell reaches r ~ 1.6r ls 5 x 10 17 cm and 
r ~ 1.6r rs 1.6 x 10 18 cm, in the thin- and in the thick- 
shell case, respectively. At this stage, ~ 90% of the energy of 
the shell has been transferred to the shocked external medium. 
Within a factor of ~ 2 in radius the blastwave fully relaxes to 
the Blandford-McKee self-similar solution. 

The different c omponents of the total energy (see e.g., 
Mimica et al.l 120071) as function of radius of the front shock 
(FS) are shown in Fig. [2] For the thin shell model, the peak 
of the thermal energy of the shell (approximately 9% of the 
total energy contained in the ejecta) traces the reverse shock 
(RS) crossing of the shell. Beyond 5 x 10 17 cm most of the 
energy that was initially in the shell has been transferred to 
the shocked external medium. The apparent "disappearance" 
of the shell at r ~ 5.4 x 10 17 cm is a numerical artifact of the 



grid re-mappingQ However, this effect is irrelevant for the dis- 
cussion of the features we are interested in, since all of them 
happen before the shell "disappears", both in non-magnetized 
and magnetized models. For the unmagnetized thick shell runs 
we see that the reverse shock dissipates more energy from the 
shell, reaching approximately 18% of the total shell energy by 
the time it crosses the shell. 

In Fig. [3] we show the relative Lorentz factor at the FS and 
the RS as function of observer time. Both shocks form immedi- 
ately after one lets the initial shell to evolve. The RS becomes 
stronger with distance, as can be seen from the increase of 
y re i with f b s . This strengthening continues until it reaches the 
rear part of the shell, where it encounters a much less dense 
medium, which results in a kink in the RS curves (marked with 
vertical black lines in Fig. [3). The peak of the emission asso- 
ciated with the RS is expected to occur at the moment the RS 
breaks out of the shell, since afterwards the density, pressure 
and velocity of the shocked medium decrease abruptly, and pre- 
cisely there it reaches its maximum strength. For the thin shell 
case this happens slightly after the burst, i.e. at observer time 
fobs ~ 1.6. For the thick shell the peak of the RS emission ap- 
pears at the end of the burst, by a time f Q b s ~ 1 -2. 

The Lorentz factor of the external medium just behind the 
FS drops as function of radius. An interesting feature is the 



2 As described bv lMimica et al.l J2007t ). grid re-mapping enables us 
to follow the evolution of a localized shell over large distance by re- 
peated re-mapping of the numerical grid. In this work the grid always 
follows the front shock, so that, once the shell has been slowed down 
by the reverse shock, it is eventually "lost" from the grid. 
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change in the slope of y^febs) at time fobs = 2.5 for the thin 
shell (at time f b s = 1.84 for the thick shell) which is a result 
of the rarefaction R2 (Fig. [TJ reaching the forward shock. For 
a short transient period the FS reduces its strength because the 
density behind it is reduced by the action of the rarefaction 
R2. The decay of y^J(fobs) is faster for Af b s ~ 0.5 after the 
rarefaction R2 catches up with the FS. Thus, we expect that the 
early afterglow emission weakens transiently faster with time. 
Later, the time dependence of y^j oc f~? expected from the 
self-similar solution is gradually established. 



4.2. Magnetized shells 

The initial phases of deceleration of strongly magnetized ejecta 
show distinct difference with respect to unmagnetized ones re- 
la ted to the magnetic pr essure of the ejecta. As we discussed 
in iGiannios et al. ( 2008 ). the jump in the physical conditions 
existent between the shell and the external medium results into 
the development of two Riemann structures, one at every ra- 
dial (rear and forward) edge of the shell. The waves into which 
each of these Riemann structures decompose are combinations 
of shocks and/or rarefactions separated by contact discontinu- 
ities. If the magnetic pressure is sufficiently large, instead of a 
typical double shock structure, a right-going shock (FS) and 
a left going rarefaction (R[) develop from the forward edge of 
the shell (Fig. |T}. Similarly, at the backward edge of the shell, 
a rarefaction (not shown in Fig. [TJ moving backwards (in the 
shell comoving frame) into the external medium develops, as 
well as a rarefaction develops penetrating the shell. 

The expansion of the shell leads to a decrease of the mag- 
netic pressure and the Lorentz factor of contact until the "tran- 
sition" radius, r tI , is reached where jqd = Jo- This distance can 
be estimated by balancing the pressure of the shocked external 
medium P - 4-"/^n e m p c 2 /3 with the magnetic pressure of the 
shell P B = B 2 /yl%7T - Eo- /87Tr 2 y 2 A (l + cr ) and solving for 
the radius 



-( 



3Ecro 



A/2 



327T«e'MpC 2 yoAo(l + CT )/ 



(7) 



From this distance onwards the shell slows down with 
respect to its initial yo. Some time after the rarefaction R[ 
has propagated backwards into the shell, a new shock devel- 
ops inside of the rarefaction fan due to the radial expansion 
of the shell (RS in Fig. [TJ. This shock sweeps backwards 
through the shell and, therefore, it is effectively a reverse shock. 
Remarkably, this shock does not immediately originate from 
the initial jump at the leading radial edge of the shell. Instead, 
it develops at the faster parts of the rarefaction fan and not di- 
rectly at the contact discontinuity separating the FS and R\ (see 
Fig. [TJ. The reason is that <x oc r 2 p in the cold, magnetized shell 
(see Appendix A), and therefore it decreases in the rarefaction 
fan. The formation of this shock can only happen when the ini- 
tial shell expands spherically, and not if the shell is assumed 
to evolve u nder conditions of pl anar (Cartesian) symmetry (as 



assumed in Mizuno et al.l 120081) . The fact that <x decreases in 



fast magnetosonic speed in the whole rarefaction fan. The de- 
crease is larger right at the head of the rarefaction where, even- 
tually, a shock formfl We therefore realize that the structure 
of the flow is much more complex than in the non-magnetized 
case, since it has developed a RS inside of the rarefaction fan 
ofRi. 

The shock is initially weak (y^ ^ 1) and remains so dur- 
ing the period in which it sweeps the whole (thin) shell (Fig.fJJ. 
When it reaches the back edge of the shell y re i 1.03, i.e., it 
remains still Newtonian. The reverse shock dissipates a neg- 
ligible amount of energy (some ~0.1% of the total energy in 
the shell). It becomes stronger in case a thick shell is consid- 
ered, reaching y re i 1 .2 by the time it reaches the rear radial 
edge of the shell (Fig. [3] lower panel). Thus, these RS sweep- 
ing a thick shell leads to a higher dissipation of energy (~ 1% 
of the total). However, the dissipated energy is still a factor of 
~ 20 lower than the in the non-magnetized thick shell model. 
In the thin shell case, the local maximum of the thermal energy 
at 1.7 x 10 17 cm (Fig. [2] upper panel) marks the RS crossing. 
Afterwards, most of the energy concentrates in the shocked ex- 
ternal medium at distance r > 4 x 10 17 cm when the evolution 
becomes very similar to that of the unmagnetized shell. 

At very early times, the Lorentz factor of the medium just 
behind the forward shock is larger than that of the shell due 
to the initial rarefaction. This leads to y re ] > 15 initially until 
the "transition" radius is reached (Fig.[3j. The fact that initially 
yrei > yo is a unique feature of magnetized ejecta (in unmag- 
netized ejecta there is always y < yo). This initial phase ap- 
pears also in the early afterglow o f the electromagnetic model 



pears also in tne early arterg 
ilLyut ikov 200(1 iGenet et all 



2007). 



The magnetization affects the (reverse) shock conditions 
and, as a result, the reverse shock crosses faster a magnetized 
shell than an unmagnetized one (see that the vertical blue lines 
appear to the left of the vertical black lines in the lower panel of 
Fig.rjJ. This fea ture has already been revealed in the study of 
(Fan et al ] |2004 . Our simulations show that the rarefaction R2 
also crosses faster (in observer time) a magnetized shell than a 
unmagnetized one. 

At a timescale a few (thick shell) or several (thin shell) 
times that of the duration of the burst the <x = and <x = 1 
models display a rather similar evolution. In this stage almost 
all the energy of the shell has been transferred to the shocked 
external medium (Fig.[2j. The two models have the same total 
energy and relax to identical asymptotic self-similar solutions 
(note the similarity between the rising parts the solid thin blue 
and black lines display in Fig. [2j. The forward shock emission 
beyond this time cannot reveal anything about the initial mag- 
netization. Note, however, that there is a rather prolonged RS 
crossing phase in the tail of the magnetized ejecta and some 
residual energy remaining in the form of Poynting flux at later 
times that may power some (energetically weak) afterglow fea- 
tures. 



the rarefaction, combined with the steepening of the rarefaction 
profile due to the spherical geometry, leads to a decrease of the 



3 We thank the referee for pointing out the fact that, without this 
decrease of the fast speed, a shock would either form immediately or 
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4.3. Dissipation by the reverse shock 

Our simulations can quantitatively answer the question of how 
much energy is dissipated when the reverse shock propagates 
into ejecta of different magnetization cr and parameter For 
practical reasons the simulations are limited to a few models0 
On the other hand, they can be used to evaluate the accuracy 



and limitations of previous (semi-)analvfical studies (|Fan et al 



12004 IZhang & Kobavashill2005l: IGiannios et alj|2008l) and use 
them as a tool to explore a larger parameter space of £ and ctq. 

The "continuous flow" models with <x = 1 and <tq — 3 
describe spherical flows of constant (as functions of radius) lu- 
minosity L, magnetization <tq and Lorentz factor yo that col- 
lide with a uniform external medium with number density n e . 
Initially the interface of the two media is set at some distance 
ro (see Section 3 for the choice of ro) and the system is let 
to evolve. With these models we can focus on the interface 
of interaction between the shell and the external medium and, 
therefore, we can track in great detail the formation and the 
strengthening of the reverse shock with time (or equivalently 
radius). This kind of set up allows us to follow the strength of 
the reverse shock for different "equivalent thickness" of shells 
in a single simulation. The idea behind the equivalent thickness 
is to measure the penetration distance from the contact discon- 
tinuity to the reverse shock and to assume a shell of the initial 
thickness Ao equal to this distance. Then we can use Eq. (HJ to 
obtain the equivalent £ of the assumed shell. The consequence 
of this is that the more the reverse shock penetrates the flow, 
the thicker equivalent shell it probes for the fixed magnetiza- 
tion <r . We define equivalent £ as (taking into account that 
E - LAq/c) 



ft q := 



31 



Ann^ipC^ 



1/6 



A- 1; V 4/3 



(8) 



one can see that thicker equivalent shell corresponds to lower 
Effectively, a run of a continuous flow model probes a line 
of constant <x in the £ - <tq plane shown in fig. S] 

A convenient measure of the strength of the reverse shock 
is the relative Lorentz factor y re i of the unshocked ejecta with 
respect to the shocked ejecta (Eq. [5]). In Fig. 21 we mark the 
locations where the y rR i becomes 1 .04, 1.1, 1.25 respectively. 
We have argued in (IGiannios et al . 2008) that for r rs <K r c , the 
magnetization of the flow cannot prevent the formation of a 
shock into the ejecta and vice versa. The curve defined by set- 
ting r rs = r c (thick black line in Fig. [U can thus be used as a 
proxy to mark the region where a reverse shock forms. 

As one can see in Fig.Hl the region of existence of revers e 
shock is larger than that estimated by IGiannios et alj (12008b . 
The reason lies on the late RS formation that has been revealed 



4 Each simulation takes between 50 and 200 thousand hours of 
computer time using between 32 and 320 processors (depending 
whether we compute thin, thick or continuous flow models) on Mare 
Nostrum (http : //www . bsc . es/plantillaA . php?cat_id=5 ) . 
The typical external storage requirements of one model vary between 
10 (thin) and 100 (thick models) gigabytes, since relatively frequent 
output of the fluid state is needed in order to obtain a satisfactory 
coverage of the fluid evolution needed for the post-processing 
calculations. 
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Fig. 4. Energy dissipation in the reverse shock in the £ - <tq 
parameter space. The circles, diamonds and asterisks mark the 
equivalent £ of the shell for which the relative Lorentz factor 
y re i of the shocked ejecta with respect to the unshocked ejecta 
becomes 1.04, 1.1 and 1.25 respectively. The dotted curves 
show the region probed by the <x = 1 and <x = 3 simula- 
tions. The solid line marks the "no reverse shock boundary" 
as estimated by (IGiannios et al ] l2008h . In the "weak/no reverse 
shock dissipation" region, the shock converts less than ~ 0.3% 
of the total energy of the shell into heat. 



by the numerical simulations. This effect was unaccounted by 
our previous analytic estimates. However, the analytic predic- 
tion that the reverse shock emission from models whose pa- 
rameters are in the region r rs ~ r c would not be observed is still 
qualitatively valid, since the dissipation from the shock is weak. 
On the solid line y re i ~ 1.1 and grows larger only for smaller 
values of For a y m \ ~ 1 . 1 , the shock converts only a fraction 
°f ~ ./bC/rei - 1)/(1 + o"o) ^ 0.1 of the en ergy of the shell into 
heat. Here ft, ~ 0.3 for o~o ~ 1 (see e.g. JZhang & Kobayashil 
20051) . Integrating the thermal energy in the shocked shell for 
the snapshot for which y re \ = 1.1, we find that it accounts only 
for ~ 0.3% and ~ 0.1% of the total energy of the shocked shell 
in the ctq = 1 and <x = 3 runs respectively. This reveals that 
the shock, though there, still dissipates very weakly "above" 
the solid line of Fig. [4] 



4.4. Reseating of the results to arbitrarily high y 

Our simulations correspond to an ultrarelativistic shell of ma- 
terial interacting with the external medium that has the qual- 
itative characteristics expected at the onset of the afterglow 
phase of the GRB ejecta. Nevertheless, they have two param- 
eters that are unrealistic with respect to what is expected in a 
GRB; namely the initial bulk Lorentz factor y an d the thick- 
ness of the shell A . GRBs are believed to originate from flows 
with yo> 100, while the thickness of the flow is likely connected 
to the observed duration of the burst through Ao ~ c?grb ~ 10 13 
cm. However, numerical reasons forced us to simulate shells 
which have yo ~ 15 and Aq ~ 10 15 - 10 16 cm. In this section 
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we demonstrate a method with which our results can reliably 
be rescaled to GRB-relevant parameters. 

4.4.1 . Motivation for the proposed rescaling 

We first focus on unmagnetized GRB ejecta and then extend 
our discussion to include magnetization. For cr = 0, the prob- 
lem of the interaction of an ultrarelativistic and thin cold shell 
with an external medium is defined by four parameters E, An, 
70 and « e . The evolution of such configuration does not depend 
on the individual parameters but on a specific combination of 
them that can be expressed as £ oc (E /n e ) 1/ ' 6 / Aq7q . The pa- 
rameter £ determines, for example, the relat ive Lorentz factor 
7,ei of the reverse shock ( Sari & Piranl 1995 ). For the ultrarel- 
ativistic flow under consideration 7 re i - ^(7sh/7o) (see Eq.O, 
where 7 s h stands for the Lorentz factor of the shocked ejecta 
when the shock reaches their rear part. This means that, for 
fixed 7 s h is a fixed fraction of 70 independently of the value 
of 7o. For example, since in our (cr ,£) = (0,0.5) model we 
have found that y s h - 0.357o - 5.2, one can predict that a non- 
magnetized shell with £ = 0.5 and arbitrary 70 » 1 is char- 
acterized by 7 s h 0.357o at the moment of the reverse-shock 
crossing. 

The idea of rescaling is to solve (numerically) for the evolu- 
tion of a shell with 70,1 and Ao,i, and then predict (without sim- 
ulating) the evolution of a second shell with 70,2 (usually larger 
than 70,1) and Ao,2 = Ao,i(yo,2/yo,i) which has the same 
£0 The distance from the central engine at which the reverse 
shock will cross the second shell is r rs> 2 = f"i- s ,i(7o,i/7o,2) 2 ^ 3 (see 
Eq. [3]>. The same relati on connects the distances r ya 
(see lsari & Piranl[l995l) where the shells "1" and "2" enter the 
deceleration phase after accumulating mass Mo/70 of their own 
initial mass Mo. This indicates that the characteristic distances 
of the shell-medium interaction for the shell "2" are shifted 
by a factor (7o,i/7o,2) 2 ^ 3 with respect to those of the shell "1". 
We postulate that the same is true throughout the evolution of 
the system. More precisely, we claim that rescaling the Lorentz 
factor from 70,1 — » 70,2 one can predict the evolution of a 
shell "2" by using that of the shell "1" providing that one also 
rescales the distance to r' — > r(Y /yo) 2 ^. 

We further extend the previous postulate by adding up the 
possibility that the shell was magnetized, i.e., we extend the 
previous claim to the case cr > 0. The evolution of a thin, 
ultrarelativistic shell with thickness 



Ao,2 = Aq, 1(70,2/70,1) 



-8/3 



(9) 



is self similar to that of a shell of the same £ and ctq and width 
A ,i- 

4.4.2. Mathematical description of the rescaling 

Here, we provide the expressions for a transformation of the 
solutions for models with the same total energy E and den- 
sity of the external medium n e . The more general transforma- 



3 For the simplicity of the discussion we fix E and n c of the two 
shells. We generalize our expressions to arbitrary E and n c in the 
Appendix B. 



tion allowing for a change of E and/or n e between two models 
is given in the Appendix B. On more mathematical terms, the 
postulated recipe for making the transformation from one so- 
lution to the other is the following. Suppose the bulk Lorentz 
factor 7i(n,fi) of the shell "1" at is known (t\ is the time in 
the rest frame of the GRB engine or laboratory frame), and we 
define the quantity / := 70,2/70,1- We further assume that the 
bulk Lorentz factors of both shells at any other time different 
from the initial one are linked by 



72(^2, h) = fy\(n,h), 
where 

h = r V3 h 

n = rps^f- 2 ' 3 + (n - r^fa)) f-«' 3 . 



(10) 



(11) 



Here r FS \ {t\) (^Vs^fe)) stands for the radius of the forward 
shock of the shell " 1 " ("2") as function of time. The other phys- 
ical quantities in the shell "2" can be derived from 72 and us- 
ing standard expressions for the forward shock conditions. The 
(gas or magnetic) pressure P in the shell and the shocked ex- 
ternal medium, the density in the shell p s heii and the density in 
the shocked external medium p ex t are given by 



Pi{r 2 ,t 2 ) = f 2 Pi(r u h) 

Psheiufo. h) = / 2 Pshell,l(ri,fl) 
Pext,2(r 2 ,t2) = /Pext,l(n,?l). 



(12) 



There are several interesting properties of the proposed recipe. 
First, the magnetization of both ejecta is the same (cro.i = cr ,2), 
not only at the initial time (which holds from our postulate), but 
also throughout the evolution, since 0"2(r2, t 2 ) °c P 3,2! Pshe\\,2 K 
Psj/Psheiu K o-\{r\,t\)- Second, the 7 re i and its time evolution 
are identical for all models (Fig. [6]upper panel), and third, the 
time evolutions of 7fs are just linearly shifted with 70 (Fig. [6] 
lower panel). 

4.4.3. Invariance of the time of observation 

An important byproduct of the transformations provided in 
Eqs. [10] - HU is that the time of observation defined in Eq. [4] 
is invariant. The only assumption we make is that the motion 
of both shells "1" and "2" is ultrarelativistic, so that the position 
of the FS can, generally, be written as 



r PS (f)=c f dT/3 FS {T)*ct-{c/2) I dry^(r). (13) 
Jo Jo 

We first demonstrate the invariance of f b s for the FS. 
Inserting Eq.[l3]into Eq.|4]we get for the shell "2" 



/■ 

Jo 



..-2, 



fobs,FS,2 = (c/2)A 2 dT 2 y/ s 2 (r 2 ) . 
' Jo 

We transform the integral as 

f' dr 2 y F 2 s 2 (t 2 ) = f-^ Pdriy^CTi), 
Jo Jo 

and insert it into Eq. [14] After transforming Ao,2 = 
/~ 8 / 3 Ao,i we finally obtain the desired result 

fobs,FS,2 = (c/2)/ 8/3 A ,ir 8/3 f ' dr irF 2 ^Ti) = f bs,FS,l - (15) 

Jo 



-2 



(14) 



Mimica, Giannios & Aloy: Deceleration of arbitrarily magnetized GRB ejecta 



I I I I I | l l l l 




r [10 ] cm 



Fig. 5. Results of a test of the rescaling with yo. Red, black 
and blue lines show the results for models with yo = 10, 15 
and 20, respectively. Full, dashed and dot-dashed lines show 
the rest-mass density, <x and the fluid Lorentz factor, respec- 
tively. Models with yo = 10 and 20 have been rescaled using 
Eqs ( [Tol l - (1121) . After rescaling of the results, the profiles of 
all models almost overlap. Only around the RS there are large 
discrepancies. The reason for them is that we have a finite time 
resolution and, therefore, we have to rescale our models using 
the closest discrete time we have to the one requested by the 
transformation expressed in Eqs. [HI 

For a point inside the shell "2", different from the FS, we 
have 

rh 

?obs,2 = (c/2)A 1 2 dT 2 y F 2 s 2O2) - cA \{r 2 - r Fsa (t 2 )) (16) 
Jo 

We see that the second term on the right hand side is also in- 
variant to the scaling, since (r 2 - ^5,2(^2)) and Ao,2 both scale 
with /~ 8 ^ 3 . This completes our proof of the invariance of f Q b s 
under the scaling relations Eqs.fTOl-fTTI 

4.4.4. Verification with test runs 

We have tested numerically the postulate stated in the previous 
section for both unmagnetized and magnetized flows and found 
that it is correct within a few percent accuracy. In the follow- 
ing, we present three numerical models which share a common 
magnetization cr = 1 and £ = 1.1. The rest of the parameters 
are: 1) y ,i = 10 and A ,i = 2.95 x 10 15 cm, 2) y , 2 = 15 and 
A ,2 = 1 x 10 15 cm and iii) 3) y , 3 = 20 and A , 3 = 2.64 x 10 14 
cm. We use the scaling relations given in Eqs. (TTOb - (fT2l to 
conform models 1 and 3 with the model 2. On Fig. [5] we show 
the density, the magnetization and the Lorentz factor after ap- 
plying the scaling laws Eqs[lO]-[T2]to the models. We show the 
scaling after the forward shock of the model 2 has reached a 
distance =s 1.75 x 10 17 cm. As we can see, Lorentz factor and 
the magnetization cr scale as expected. The rest-mass density 
within the shell and the pressure (not shown here) are also fol- 
lowing Eqs. Q~2] In Fig. [5] we note that there is a factor of 2 
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Fig. 6. Similar to Fig. [3] but for the models used in the test of 
the rescaling hypothesis. Black, blue and red lines in the upper 
(lower) panel show the relative Lorentz factor of the fluid at the 
FS (RS) as a function of the normalized time of observation for 
models with y () = 10, 15 and 20, respectively. As expected, the 
relative Lorentz factor at the front shock scales with yo while 
at the reverse shock they coincide. The dashed line in the upper 
panel shows the results of a model obtained by rescaling to 
yo = 100 the model with yo = 15. 



difference between the rescaled results around the RS. The rea- 
son for this discrepancy arises from the finite time resolution 
of our models. In order to rescale them we have to use the clos- 
est discrete time we have available to the one requested by the 
transformation expressed in Eqs.fTTI 

4.5. Extrapolation to y > 100 and "bolometric" light 
curves 

The MHD calculations presented here do not suffice to calcu- 
late the expected emission in detail. Such a calculation of the 
afterglow emission in different observed bands needs additional 
assumptions related to the, poorly known, shock microphysics. 
These include, for example, the energy distribution of the ac- 
celerated electrons and the generation of magnetic fields in the 
shock front. Furthermore, one needs to include radiative mech- 
anisms such as synchrotron and inverse Compton and allow 
for adiabatic losses. This detailed calculation falls beyond the 
scope of this work. Instead, we apply a simple method to get 
a rough estimate of the bolometric emission expected from the 
different models. In order to make predictions about the prop- 
erties of the afterglow light curves associated to our numerical 
models, we extrapolate them to the conventionally accepted pa- 
rameter regime where GRBs shall take place. To do this, we 
apply the rescaling described in the previous Section. 

We assume that in both reverse and forward shocks a frac- 
tion e e « 0.1 of the dissipated energy goes into high-energy 
electrons, and that the electrons are in the fast cooling regime 
(as is usually the case during the initial afterglow phases). Then 
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the resulting total emission is given by the rate at which the 

shocks heat the external medium and th e shell. 

From the shock conditions at the FS ( Sari & Piranll995 ) we 
find that the thermal energy (in the laboratory frame) dissipated 
by the FS when it moves from r to r + Ar is 



t [seconds] 



A£ t h(r, t) — \6nr 2 Ary 2 n e m p c 2 



(17) 



where we assume that Ar <K r, and that y is approximately 
constant between r and r + Ar. We assume the luminosity of the 
FS to be a fraction of the thermal energy dissipated due to the 
heating of the external medium by the shell in a time interval 
At = Af bsAo/c, where Af b s is the difference of normalized 
observational time (Eq. H} between the moments in which the 
FS moves from r to r + Ar. According to this definition, the 
bolometric luminosity for the front shock reads 



Lps := 16nn e m p c 2 e e y 2 r 2 



Ar 
Ar 



(18) 



For the reverse shock we estimate its luminosity assum- 
ing that a fraction e e of the increase of thermal energy in the 
shocked shell, as it moves from r to r + Ar, is radiated by the 
RS. Thus, we define 



Lrs 1= (Af) max 



0,Q f 

Jsr 



4-y 2 p - I 4y 2 /7 

shell r+Ar v/ shell r 



, (19) 



where the f. „ denotes the integral over the volume of the 

Jshell r ° 

shell when the FS is at r. The integrand 4-y 2 p can easily be 
derived from the e quations of RMHD assuming an adiabatic 
index 4/3 (see e.g. jMimica et al ] l2007l) . 

The luminosities Lps and Lrs can be normalized to the ini- 
tial shell luminosity defined as 



L Q := 47rro7op shell c 3 . 



(20) 



In this way we obtain the dimensionless luminosities /ps := 
L FS /L and l RS := Lrs/Lq- 

We have checked that the normalized and conveniently 
scaled light curves for the test models from Sec. |4.4.4| coincide 
to within a few percent deviation. This means that we can use 
such a normalized light curves to predict what would be the 
bolometric luminosity observed from a shell which has, e.g., 
yo = 100, and the same value £ = 1.1 and 0.5 as our thin and 
thick shell models, respectively. The light curves computed for 
the reference models with yo = 15 and scaled to ejecta with 
Lorentz factor yo = 100 (thin shell; Fig. [7} and yo = 300 (thick 
model shell; Fig. [8} show a fundamental difference between 
thin and thick magnetized shells. The luminosity of the RS of 
magnetized thin shells is much smaller than the luminosity of 
the corresponding RS in the hydrodynamic case (note that in 
Fig.[7]the light curve of the magnetized RS does not even show 
up at the scale we are considering). For thick shells, the lumi- 
nosity of the RS shock is about 10 times smaller than that of 
the corresponding thick shell with <x = 0. Hence, the detec- 
tion of the RS will be, in general, much more difficult if the 
shell ejecta is magnetized than if it is unmagnetized. Indeed, if 
the magnetized ejecta is thin, it is very likely that the RS is not 
detected at all. 




Fig. 7. Light curves for non-magnetized (black) and magne- 
tized (blue line) thin shell models after scaling the yo = 15 
model to yo = 100. The luminosity of the FS is shown in thick 
lines, while the black dashed line shows the luminosity of the 
reverse shock of the un-magnetized shell. The RS of the mag- 
netized shell has a very weak dissipation and its light curve 
is not shown on the scale of the plot. Left and bottom axes 
show the dimensionless luminosity (Zps and irs) and the time 
of observation, respectively. The displayed luminosity profile 
can be applied to any other shell with a different yo provided 
that £ = 1.1 (see text). Right and top axes show respectively 
the luminosity and the normalized time of observation for the 
particular case yo = 100. 
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Fig. 8. Same as Fig. but for a thick shell model = 0.5) 
at yo = 300. The reverse shock emission from the magnetized 
thick shell (dashed blue line) is much stronger than in the mag- 
netized thin shell. 



5. Discussion/Comparison with previous work 

GRB outflows may be launched by strong fields resulting in 
a Poynting-flux-dominated wind. In this case the flow can re- 
main strongly magnetized throughout the acceleration, colli- 
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mation and GRB emission phases all the way to the onset of 
the afterglow. This is in sharp contrast to the expectations from 
a flow initially dominated by thermal energy (fireball) in which 
magnetic fields are dynamically unimportant. Magnetization is 
expected to affect the initial phase of interaction of the GRB 
ejecta with the external medium. Since early afterglow obser- 
vations are now possible for many bursts, it is becomes topical 
to study the effect of magnetization of the ejecta in the early 
afterglow in more detail. 

Here we perform ultra-high resolution ID RMHD simu- 
lations of the interaction of a radially expanding magnetized 
ejecta with the interstellar medium. We study the complete evo- 
lution of an ultrarelativistic, initially uniform ejecta shell all the 
way to a the self-similar asymptotic phase. We explicitly show 
which is the required resolution of numerical simulations to re- 
solve appropriately all the discontinuities in the flow, and to be 
free of numerical artifacts. The main findings of the work are 
the following: 

1. The complete evolution can be characterized in terms of 
two parameters, namely, the Sari-Piran parameter £ and the 
magnetization cr . Since both parameters are defined from 
combinations of basic physical properties of the ejecta (e.g., 
yo, E, Ao, etc.), a single point in the {£, cr ( )}-plane can be 
used to probe a variety of equivalent combinations of basic 
physical parameters. A very useful byproduct of this de- 
generacy is that one can produce numerical models where 
the shell Lorentz factor is yo ~ tens and rescale the results 
to arbitrarily large yo. Certainly, the numerical difficulty of 
simulations with moderate Lorentz factor (namely, yo < 30) 
is smaller than those with a large one (see App. IA.lT ). The 
later type of simulations are prone to a number of numerical 
artifacts that hinder an appropriate comprehension of the 
physics we are dealing with. The method to rescale from 
our reference numerical models to the conditions expected 
to occur in GRB afterglows is described in Sect. 14.41 and in 
App. B. 

2. The amounts of energy dissipated in the reverse shock de- 
pend strongly on the magnetization of the flow. The reverse 
shock is very weak or absent for ejecta characterized by 
£ > 0.5 (thin shells) and o~o > 1 . The emission from the re- 
verse shock is strongly suppressed for <t () > 1 except for 
high yo > 1000 flows (or equivalently low £ flows). More 
moderate values of the magnetization cr ~ 0.1 can lead to 
powerful emission, in excess to that expected from cr = 
ejecta, since there is both a strong shock and a strong mag- 
netic field for efficient synchrotron emission. The last state- 
ment assumes that there is particle acceleration taking place 
in magnetized shocks with efficiency similar to that of un- 
magnetized shocks. 

3 . Magnetically dominated (cr > 1 ) ej ecta are characterized by 
an initial rarefaction originated at the leading radial edge of 
the shell that results in a Lorentz factor of the front shock 
y > yo- The energy that is dissipated by the forward shock 
and the expected emission at the onset of the afterglow is 
much higher for cr > 1 ejecta with respect to weakly mag- 
netized ones. "Bolometric" light curves are presented in 
Sec. 1431 



The magnetic energy of the shell is increased due to shock 
compression during the reverse shock crossing in a gree- 
ment to the findings of Zhang & KobayashU ( 2005 ). We 
have found that the bulk of the magnetic energy is trans- 
ferred into the shocked external medium on a rather sort 
timescale (a few to several times the burst duration) for the 
£ = 1.1 and £ = 0.5 models we have simulated. Several 
light crossing times of the initial ejecta width suffice for the 
evolution of ejecta of cr ~ 1 to become very similar to that 
of the cr = simulation. At this stage almost all the en- 
ergy has been transferred to the shocked external medium 
and the forward shock emission is practically the same in- 
dependently of the initial magnetization of the flow. One 
should look to the onset of the afterglow to identify effects 
connected to the magnetization. 



Comparing with previous work, Zhang & KobayashU 



(2005) significantly overestimated the reverse-shock emission 
for (To ~ 1, mainly because of the assumed higher value of 
the relative Lorentz factor y rR i (and conse quently dissipation) 
in the reverse shock. The iGiannios et al.l analytical curve on 
the {£, (To)-plane provides a good approximation of when there 
is substantial reverse shock dissipation. The conclusion of that 
work that the observed paucity of optical flashes -signatures of 
reverse shock emission predicted by the fireball model- may be 
understood by <r > 1 ejecta is verified by our simulations. 

Modeling of the emission associated with the forward and 
reverse shock can be used to compare the magnetization of 
the s hocked external m edium to that of the shocked ejecta 



(e.g, Zhang et al. 2003b. This method has been applied in a 



number of bursts dFan et aL 20021: Kumar & Panaitescu 2003 



Mundell et all l2007bl: iGomboc et al.ll2008l) . We caution here 
that such approach considers hydrodynamical shock conditions 
and is therefore not reliable when the magnetization of the 
ejecta is large enough to alter the dynamics of the decelera- 
tion. In the latter case a full MHD calculation (such as the one 
presented here) is needed for fruitful comparison to observa- 
tions. 

There are aspects of the problem of interaction of magne- 
tized ejecta with the external medium that have not been settled 
by this work. Although we solved for the dynamical evolution 
of the ejecta, the strength of the shocks as function of time 
and computed approximate "bolometric" light curves, we did 
not calculate detailed light curves in different observed bands. 
For this calculation additional assumptions related to the shock 
microphysics and inclusion of radiative processes such as syn- 
chrotron and inverse Compton emission are needed. 

Furthermore, we have not explored the high £ (or 
Newtonian reverse shock) regime. In this regime, shell spread- 
ing because of the presence of slower and faster parts within 
the shell has to be considered. For £ > 1, the onset of the af- 
terglow takes place on a later observer time and can be used 
to infer physical quantities such as the Lorentz f actor of the 
flow y dSari & Pirarilll999t iMundell et al.ll2007ah . The slope 
of the initial rising part and the peak of the curve depend on 
the external medium density profile and probably the magneti- 
zation of the ejecta. These features are worth to be investigated 
in more detail. Finally since the initial interaction lasts longer 
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as seen by the observer, substantial magnetic energy remains in 
the shell. It is possible that dissipation of this magnetic energy 
is localized active r egions results in late-time flares as proposed 



Giannios! (12006). 
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Appendix A: Numerical method 

We solve equations of RMHD in ID spherical geometry as- 
suming that the fluid moves only in the radial direction. The 
magnetic field is purely toroidal magnetic, i.e., the magnetic 
field which is perpendicular to the direction of propagation of 
the fluid. The system of RMHD equations is (with the speed of 
light set to be c — 1) 



dt r 2 dr\ I 



(A.l) 



where the vector of unknown or conserved variables are 

U = (jjy, ph*y 2 v, ph*y 2 - p* — py, fi) . (A. 2) 

The fluxes in Eq. lA.ll are 

F = (pyv,ph*y 2 v 2 + p* , [ph'y 2 - py] v, vBj , (A. 3) 
and the source terms read 



(A.4) 



Here p, p, y and B are the fluid rest mass density, pressure, 
Lorentz factor and magnetic field in the frame of the central 
engine or laboratory frame. The magnetic field B is measured 
in Gaussian units. The total pressure is p* :- p + B 2 /2y 2 , and, 
the specific enthalpy is h* := 1 + yp/(y - l)p + B 2 /py 2 . In our 
models the fluid is assumed to be an ideal gas with the adiabatic 
index y = 4/3. We note that we express the components of all 
three-vectors in the physical, i.e., orthonormal basis. 
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Romero et alj ( 12005b discuss the solution to the Riemann 
problem in case in which the magnetic field is perpendicular 
to the fluid velocity and in Cartesian geometry. They show that 
the ratio B/(yp) is constant everywhere except across contact 
discontinuity. An analogous expression in spherical geometry, 
cr oc rp can be derived assuming a cold magnetized fluid. In 
this case the system of equations I A. 1 1 reduces to three equa- 
tions. From the continuity and the induction equation one can 
easily derive the desired relation. 

We use t he relativistic magnetoh ydrodynamic code 
MRGENESIS dMimica et al.l 120051 120071) . a high -resolution 
shock capturing sche me 
1999: 



based on GENESIS dAlov et al 
Leismann et al.l 12005b In our code the fluid is dis- 
cretized in spherical shel ls (zones). We use the PPM 



dColella & Woodward! 119841) scheme for the spatial interpo- 
lation of variables wit hin numercial zones, and a HLLC 
dMignone & Bodoll2006l) approximate Riemann solver to com- 
pute numerical fluxes accross zone boundaries. The time inte- 
gration is performed using a third-order Runge-Kutta method. 

A .1. Numerical resolution 

For the simulation results to be as free as possible of numeri- 
cal artifacts a large enough resolution is needed. Of particular 
concern is the minimum number of iterations from the start of 
the simulation which are necessary to resolve the initial evo- 
lution of the discontinuity that separates the shell from the ex- 
ternal medium (that forms a Riemann problem). The knowl- 
edge of this information is requiered since almost all RMHD 
codes based on approximated Riemann solvers develop ini- 
tial transient spureous behaviours at the location of the orig- 
inal discontinuity. These spureous behaviours relax with time 
to the correct physical solution. Therefore, since the problem 
is self-similar in Cartesian coordinates, almost independently 
of the initial resolution, our numerical code recovers correctly 
(i.e., within the accuracy of our method) the physical solu- 
tion. However, in spherical symmetry the problem is not strictly 
self-similar. Thus, the break up of an initial discontinuity may 
yield to the formation of additional discontinuities (inside of 
the Riemann fan but not directly emerging from the contact 
discontinuity) at later times. This is precisely what happens in 
the rarefaction R i in our magnetized models (see § I4.2I ). where 
the RS forms. If the formation of the shock happens very close 
to the location of the contact discontinuity, the initial transient 
artifacts in numerical simulations may pollute the formation of 
the RS and yield to a wrong numerical solution, where, e.g., the 
RS does not form. The way to diminish the hampering effect of 
such initial transients is to increase the numerical resolution 
around the initial discontinuity. 

To eliminate the effects of the spherical geometry on the 
solution of the Riemann problem, and in view of the fact that 
our initial shells start at distances Rq > 10 16 cm, where the local 
effects of the spherical geometry are practically negligible, in 
this section we study both exact and numerical solutions of the 
following Riemann problem in planar coordinates: 

- left state: p = 1, p = 10~ 2 , yo = 15 and B = 15 , 

- right state: p = 10~ 4 , p = 10~ 6 , y = 1 and B = . 




Fig. A.l. Time t (upper panel) and the number of iterations Ai ter 
(lower panel) needed to resolve the Riemann problem in planar 
coordinates as a function of the spatial discretization Ax. 
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Fig. A.2. Time needed to resolve the Riemann problem in pla- 
nar coordinates as a function of the initial Lorentz factor Fo for 
Ax = 1.5625 x 10~ 5 . 



The solution to this Riemann problem is self-similar, and con- 
sists of a right-going shock wave separated by a contant discon- 
tinuity from a left-going rarefaction wave. We are interested in 
the time t and the number of iterations Ai ter it takes a numerical 
code to obtain a correct Lorentz factor (to an accuracy of less 
than a percent) of the contact discontinuity (yen = 25.56 for 
this particular problem). 

Figure I A. 1 1 shows the results of the test for seven differ- 
ent simulations whose zone sizes have values Ax = 10~ 3 , 5 x 
10~ 4 , 2.5 x 10~ 4 , 1.25 x 10~ 4 ,6.25 x 10~ 5 , 3.125 x 10~ 5 and 
1.5625 x 10~ 5 . By fitting to the data points in the plot, we 
find r oc Ax 1 003±0002 . This linear dependence can be seen i 
the lower panel of Fig. IA. H where we see that A; te r is roughly 
independent of the resolution (Fig. IA.ll lower panel). We find 
that t oc Ax A; ter . 

Of equal importance is the dependence of r and Ai ter on 
the initial Lorentz factor. We modify the Riemann problem by 
changing yo of the left state and the magnetic field, such that 
the ratio B/yo - const. For this test we use the finest resolu- 
tion Ax = 1.5625 x 10 for all models. Results are shown on 
Fig. IA.2I We find r oc yi-08±o.02 Due to t h e constant Ax this 
means that Ai ter also depends almost linearly on yo. Combining 
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results shown in Fig. lA.ll with those of Fig. |A.2| and assuming 
linear dependences of t on Ax and yo we conclude that the re- 
quired spatial discretization for Riemann problems similar to 
those discussed in this paper follows the relation 

Ax oc tjq 1 . (A.5) 

It has to be pointed out that the constants which are implicit 
in Eq. I A. 5 1 depend on the initial density and pressure ratio, as 
well as the magnetization of the fluid. The result expressed by 
Eq. |A.5| can also be interpreted in the following way: for a fixed 
numerical resolution, if the Lorentz factor of the problem to be 
solved grows, the time needed to relax any initial numerical 
pathology also grows. 

We make use of Eq. IA.5l to determine the maximum Ax (or, 
conversely, the minimum resolution) needed to run our simu- 
lations such that t is much smaller than any of the characteris- 
tic hydrodynamic time scales of our models. Particularly, t we 
warrant that t is smaller than the time needed to form the RS 
in the rarefaction fan Ri . 

A reduced density jump between the shell pl and the exter- 
nal medium pr reduces drastically the numerical complexity of 
the break up of the Riemann problem. Our choice of the den- 
sity jump Pl/pr ~ 10 4 tries to reach the large density contrast 
expected in the conditions found in GRB afterglows (although 
it is still smaller than what an optimal modeling demands). We 
point ou t that a much reduc ed value Pl/Pr ~ 10 2 (as consid- 



ered by iMizuno et al.l 12008) could be too small in regard of 



conditions to be met in this astrophysical context. 



